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ABSTRACT 



This thesis compares the outputs of the ionospheric 
propagation prediction model ADVANCED PROPHET, 4.3 to 
measurements of propagation mode delay for a High Frequency 
communications link between Monterey and San Diego, 
California. 

Mode delay variations throughout the day are presented for 
experimental data and PROPHET predictions. A margin of error 
of less than 0.5 msec was considered acceptable and the number 
of acceptable predictions per day was generated. Acceptable 
predicted data collected over the test period was analyzed to 
establish which hours of the day PROPHET accurately predicts 
propagation mode delay, independent of frequency, date and 
power levels. During the first six hours of the day PROPHET 
data tracks experimental data for mode delay change patterns. 
On a daily basis, predictions are best between 1400 and 1700 
GMT (0600 and 0900 local time) , although patterns could not be 
established for other hours of the day. 

Predicted mode delay percentage distributions show greater 
correlation during the first three hours of the morning and at 
sunset . 
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I . INTRODUCTION 



The concept of an electrical conducting layer in the 
atmosphere dates back to the mid- 19th century. Shortly after 
the invention of radio, when communications distances began to 
increase, the relation between propagated waves and this 
electrical layer, now known as the ionosphere was observed. 
Since then, many studies have been conducted on the 
characteristics of the ionosphere and its effects on radio 
communications . 

The principal use of the ionosphere is in the High 
Frequency (HF) band. Since communication in the HF band 

depends on ionospheric wave reflection, which is continuously 
changing due to daily, seasonal, and solar conditions, the 
effect of the changing ionosphere on radio communications must 
be considered. Ionospheric communication, when reliable, has 
several advantages over other methods. The low cost, 
simplicity, and portability of equipment are advantages over 
other communications methods that may be more reliable. 
Communications via the ionosphere is wide spread among CIS and 
third-world countries and is being given consideration as a 
viable alternative by other users. 

It is impossible to maintain reliable communications 
without considerating the changing environment in which it 
occurs. Hence, it is important to have knowledge of 
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ionospheric characteristics and the elements that affect it 
such as time of day, season, location, and solar activity [1] . 

Computer models permit the analysis of HF wave propagation 
parameters and can predict propagation channels that will 
support communications. 

A project named Polar Equatorial Near- Vertical Incidence 
Skywave Exp eriment (PENEX) which acquires calibrated field 
strength measurements from ionospheric radio transmissions is 
creating a data base for comparison with different ionospheric 
propagation prediction models. In addition, as part of the 
PENEX project, a communication link has been established for 
the purpose of equipment testing, calibration and experimental 
methodology . 

After this communication link between Monterey, CA. and 
San Diego, CA. is tested, project PENEX will be expanded to 
include a transmitter site in Alaska, and receiver sites in 
Alaska, California, Pennsylvania and Washington. 

This thesis uses one of the computer models, ADVANCED 
PROPHET 4.3 to predict HF skywave propagation for comparison 
and correlation to the experimental data obtained for mode 
delays in the communication link between Monterey and San 
Diego . 

It is also important to note that due to limited data 
availability, the research analysis of this thesis is somewhat 
restricted. 



2 



II. THEORETICAL BACKGROUND 



A. THE IONOSPHERE 



The 


ionosphere 


is 


the 


region of 


the atmosphere 


which 


extends 


between 50 


to 


600 


km where 


sufficient ionization 


occurs 


in layers 


to 


affect radio 


transmissions . 


The 



ionosphere is an electrically neutral region, since if it had 
a given excess charge, then electrical forces would prevent 
the formation of stable layers [2] . 

1. Photo- ionization 

Photo- ionization occurs when the sun produces extreme 
ultraviolet (EUV) light waves which in turn detach electrons 
from neutral atoms in the Earth's atmosphere. When this takes 
place, the atoms become positive ions as illustrated in Figure 
1 . 

Since the atmosphere is bombarded by EUV of different 
frequencies, different layers are formed at different 
altitudes; higher frequency EUV waves penetrate the atmosphere 
deeper, producing layers at lower altitudes while lower 
frequency EUV penetrates less, producing higher altitude 
layers [3] . 

2 . Recombination 

When collisions are produced between positive ions and 
free electrons, neutral atoms are formed. Recombination can be 
seen as the photo- ionization process in reverse. 
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Figure 1. Photo-ionization of a neutral atom, A, by extreme 
ultraviolet light from the sun, yielding a 
positively charged ion, A*, and a free electron, e. 
[Ref.l, p. 18] 
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Between dawn and dusk, photo- ionization exceeds 
recombination and maximum density layers are formed affecting 
radio communications in greater percentages than during the 
hours after dusk when the effect is reversed. 

3 . Ionospheric Layers 

The ionosphere divides into three mayor layers, D, E , 
and F which is sub-divided into two others, FI and F2 . 

The lowest layer is the D- layer which spans from 50 to 
90 km. In the D- layer low frequency (LF) waves are refracted 
while Medium Frequency (MF) and High Frequency (HF) waves are 
absorbed. There are two types of absorption that affect the 
HF band: deviative and nondeviative . Deviative absorption 
occurs when the refractive index approaches zero, generally at 
the apex of the trajectory, while nondeviative absorption 
takes place when the refractive index approaches one and the 
product between electron density and electron collision 
frequency is high [2] . Another aspect of the D- layer is that 
it disappears after sunset due to rapid recombination. 

The E-layer is located between 90 and 130 km and is also 
known as the "Kenelly-Heaviside layer”. This layer is 
characterized by strong diurnal variations. Recombination is 
rapid after sunset and the layer diminishes as the night 
advances. The Sporadic E- layer also occurs near 120 Km and 
exists throughout the day at lower (equatorial) latitudes and 
during late night for higher (polar) latitudes. The E-layer 
reflects signals up to 20 MHz providing communications at 
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distances up to 2500 km per reflection. 

The F- layer is formed between 150 and 600 km. During 
the daytime this layer divides into two layers, FI and F2 , 
located between 150 and 220 km and above 225 km, respectively. 
Since these layers are at the highest altitudes, photo- 
ionization occurs at the highest rate. At noon time in the 
northern latitude winter, these layers are the closest to the 
sun and the ionization rates are maximum. Because the 
atmospheric density diminishes with altitude, the F- region 
recombination process is relatively slow and at night the 
F- layer is always present. Figure 2 shows the different 
layers of the ionosphere. 

4. Variations of the Ionosphere 

Since the ionosphere is created by the effect of solar 
radiation upon the atmosphere, the ionosphere will also vary 
with the time of the day, season, location on the earth, solar 
activity, and other similar factors. 

These variations that affect the ionosphere will also 
affect HF communications. The principal variations of the 
ionosphere are: 

* Diurnal (throughout the day) 

* Seasonal (throughout the year) 

* Location (geographic and geomagnetic) 

* Solar Activity (Solar cycle and disturbances) 

* Height (different layers ) [1] 

Figure 3 illustrates how diurnal variations affect the 
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Figure 2. Sketch of the various regions of the ionosphere 

as it appears during the day and during the night. 
[Ref. 1, p. 17] 
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Figure 3. Diurnal variation of electron density and the 

formation of ionospheric layers. [Ref. 6, p. 21] 



8 



formation of ionospheric layers. 

5. Ray Tracing 

In order to determine the amplitude, phase, 

polarization, flight time, etc. of the transmitted wave, it is 
necessary to determine the ray path between the transmitter 
and the receiver. Determining the ray path seems to be a 
simple task but in practice it is not. Factors affecting the 
ray path are the orientation and intensity of the earth's 
magnetic fields, the collision frequency of electrons and 
electron density variations within the layers. Thus, 
calculations are much more complex for the real-world than for 
the assumption that the ionosphere is a homogeneous layer. 

6 . Ray Paths in the Ionosphere 

The path followed by a wave reflected from the 
ionosphere depends on the incidence angle at the ionosphere. 
After entering the ionosphere, a wave can take different 
paths. It could reach the receiver via one hop, i.e., 
bouncing off the ionosphere once, or through multiple hops. 
Reflection between the different layers in the ionosphere as 
well as any combination of the above cases is also a 
possibility. Figure 4 shows examples of different paths for 
a wave using the ionosphere. 
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Figure 4. Raypaths for a fixed frequency. [Ref. 



III. THE MONTEREY-SAN DIEGO COMMUNICATION LINK 



Before the polar region HF ionospheric communication link 
for PENEX is installed, a test-bed link was established to 
verify the operation of the equipment and software. For this 
link, Monterey was designated as the transmitter site with San 
Diego as the receiver site. The geographical latitude and 
longitude used in the scenario files of ADVANCED PROPHET 4.3 
are 36°37' N latitude and 121°55' W longitude for Monterey and 
32°43' N latitude and 117°09' W longitude for San Diego. The 
transmit power for the tests is 2.5, 75 and 100 watts and the 
frequencies are 5.604 and 11.004 MHz. The methodology of 
experimentation is described in detail in Chapter IV. 

A. LINK ANTENNAS 

The description of the antenna characteristics for the 
transmitter and receiver sites are essential for modeling the 
skywave propagation and for the link power budgets. 

1. Transmitting Antenna 

The transmitting antenna utilized is a multiband dipole 
which consists of a group of three half -wavelength dipoles 
connected to a common transmission line. The multiband dipole 
is designed to operate on 5 . 6 , 11 and 16.8 MHz respectively. 
For data collection purposes the antenna was operated at 5.604 
and 11.004 MHz since 16 MHz propagation was very unlikely for 
the test link. The data collected for the 11.004 MHz 
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frequency range was obtained for only one day, during the time 
of this investigation. 

Dipole lengths were 5.6, 11, and 16.8 MHz for 24.69, 

13.05, and 8.53 meters, respectively. Since these dipoles are 
center- fed, each has two elements half the dimensions given 
above. Figure 5 shows the dipole antenna geometry details. 
The separation between the dipoles is 0.3048 meters and the 
upper and lower elements are at an elevation angle of 11.21° 
and -11.21° respectively. The dipoles are fed by a 50 ohm 
impedance coaxial cable connected to a 1:1 balun transformer. 

2 . Receiving Antenna 

The receiver antenna in San Diego is a dipole, 
containing three hal f - wavelength dipoles on one structure, 
drooping from the center at -34°. It is designed to operate 
at the same frequencies as the transmitting antenna. 

The half lengths of each element are 12.73 m. , 6.48 
m., and 4.24 meters for frequencies 5.6, 11, and 16.8 MHz 

respectively. The separation between elements is 12.7 cm and 
the feed line is a 50 ohm coaxial cable. Figure 6 shows the 
details of the receiver antenna. 

3 . Antenna Simulation 

The transmit antenna designs were modeled by a Naval 
Postgraduate School student using the Numerical 
Electromagnetics Code (NEC) which in turn provided patterns 
for ADVANCED PROPHET 4.3. Figures 7 and 8 show radiation 
patterns from NEC. In ADVANCED PROPHET 4.3, these patterns 
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Plgura 5. Montaray Multiband Dipola Antanna 
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Figure 6. San Diego Drooping Dipole Antenna 
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are required only for elevation angles of 6, 20, 40, 50, 70, 
and 90 degrees. 

Since radiation pattern definition in ADVANCED PROPHET 
4.3 is for the entire HF band and for only integer values of 
frequency, the pattern for 5.6 MHz was used for 2, 3, 4, 5, 6, 
7, and 8 MHz. Likewise, the 11 MHz pattern was used for 9, 
10,11,12 and 13 MHz, while the 16.8 MHz radiation pattern was 
used for 14, 15, 16, 17, 18, 19, 20, 25, and 30 MHz. Figures 
9 and 10 show the radiation patterns and antenna 
characteristics used by ADVANCED PROPHET 4.3. 
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ELEVATION PATTERN AT 0 
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Figure 7. Elevation radiation patterns obtained by NEC for 
the Multiband Dipole Antenna 
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the Drooping Dipole Antenna 
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Figure 9. Elevation radiation patterns used by PROPHET for 
the Multiband Dipole Antenna 
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Figure 10. Elevation radiation patterns used by PROPHET for 
the Drooping Dipole Antenna 
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IV. COMPARISON OF EXPERIMENTAL DATA WITH PROPHET PREDICTIONS 



A. EXPERIMENTAL DATA 

The following data are collected via PENEX equipment and 
software at the San Diego receiver site. The PENEX system 
uses direct sequence spread spectrum modulation (DSSS) , 
specifically phase shift keying (PSK) , to transmit a reference 
maximum length sequence. The system measures mode time delay 
and signal strength using a quadrature - sampled, Fast Fourier 
Transform (FFT) matched- filter signal processor. The time 
delay measurements have a 25 millisecond (msec) resolution. 
This system provides a precise time delay measurement as well 
as consistent signal strength measurements. 

Figure 11 shows an example of correlation vector plots 
resulting from cross correlation of the reference sequence 
with the digitized HF data. The upper plot is the raw vector 
(linear form) , of approximately 10 msec of data. Since the 
digitizing rate is 39.0625 MHz, each bin is 25.6 microseconds. 
In order to determine the mode delay, delta t, and the value 
of the correlation peak, this vector is scanned to locate the 
peaks, record the corresponding correlation peak values, and 
find the distances between adjacent peaks (in terms of the 
intervening bin count) . The number of intervening bins is 
then multiplied by the 25.6 microsecond duration of each bin, 
yielding the mode delay, delta t. Next, the value of the 
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Figure 11. Correlation vector plot that results from the 

cross correlation of the reference sequence with 
the digitized HF data. 
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correlation peak is mapped, through a calibration curve, to 
received voltage and later, through the CCIR field strength 
equation, to field strength. It is important to note that the 
dB value obtained for all three modes in Figure 11 is not the 
field strength but rather a processing gain measurement 

(20*LOG{max correlation/noise floor}). 

The information produced by the PENEX team consists of mode 
delay vs. real time graphs for five different days, as shown 
on Table 1; three different transmission power levels and two 
frequencies were used. 



TABLE 1. DATE OF TRANSMISSION, FREQUENCY AND POWER USED 
BY THE MONTEREY- SAN DIEGO COMMUNICATION LINK. 



Date 


Frequency (MHz) 


Power (watts) 


28 Aug 92 


5 . 604 


2.50 


02 Oct 92 


5 . 604 


75.00 


03 Oct 92 


5 . 604 


75.00 


04 Oct 92 


5 . 604 


75.00 


07 Oct 92 


5 . 604 


75 . 00 


| 02 Oct 92 


11 . 004 


100 . 00 



The plots shown in Figures 12 through 17 provide time delay 
for 24 hours of each day. For the purpose of comparison, the 
data considered were the data corresponding to each hour 
interval. All mode delays were stored on file for later 
comparisons. All times are in GMT, (the Monterey-San Diego 
Link has a -8 hour offset) . Table 2 shows the time 
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Figure 12. Measured Mode Delay for 28 Aug 92 
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Figure 13. Measured Mode Delay for 2 Oct 92. 
(5.604 MHz - 75 watts) 
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Figure 14. Measured Mode Delay for 3 Oct 92. 
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Figure 15. Measured Mode Delay for 4 Oct 92. 
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Figure 16. Measured Mode Delay for 7 Oct 92. 
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Figure 17. Measured Mode Delay for 2 Oct 92. 
(11.004 MHz - 100 watts) 
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correspondence between GMT and link local time. 



TABLE 2 . TIME CORRESPONDENCE BETWEEN 



GMT AND LOCAL TIME 



GMT 


Local Time 


GMT 


Local Time 


00 : 00 


16 : 00 


12:00 


04 : 00 


01:00 


11:00 


13 : 00 


■ 

0* : 00 


02 : 00 


18:00 


14 : 00 


06 : 00 


03:00 


14:00 


15 : 00 


07:00 


04 : 00 


20:00 


16 : 00 


08:00 


05:00 


2§ : 00 


17:00 


04:00 


o 

o 

o 


2§ : 00 


18:00 


• 6:00 


0« : 00 


23 : 00 


13:00 


11:00 


08 : 00 


116 : 00 


211: 00 


12:00 


1 09:00 


01:00 


21:00 


13:00 


• 0:00 


02:00 


22 : 00 


• 4 : 00 


11 : 00 


03 : 00 


23:00 


15:00 



B. ADVANCED PROPHET 4.3 RESULTS 

1. Input Scenario Parameters for PROPHET 

In order to have a controlled experiment and a basis for 
comparison between experimental and computer predicted data, 
an environment similar to the experimental scenario was 
created. The first step in using ADVANCED PROPHET 4.3 is to 
define the transmitter and receiver antenna characteristics 
and receiving site noise. As mentioned in Chapter III the 
radiation patterns for both antennas were obtained via NEC 
simulation. Man-made noise was defined as QM-type (quasi- 
minimum) . The 10.7 cm. solar radio flux as well as Kp (the 
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planetary magnetic index) were provided by Mr. Wayne Bratt of 
the Naval Research and Development Center (NRAD) San Diego, 
shown in Table 3. Since Kp is an average of the variation of 



TABLE 3. MAGNETIC INDEX AND 10.7 cm. SOLAR RADIO FLUX 
USED IN PROPHET SIMULATIONS. 




the magnetic field on the earth over a three-hour period, 
eight (8) constants for each day were divided in three hours 
intervals. RAYTRACE, a program that is part of ADVANCED 
PROPHET 4.3, was used to obtain the different time delays for 
the Monterey-San Diego communications link. This program 
calculates ray paths that propagate between the transmitter 
and receiver antennas for a given frequency, time and 
ionospheric parameters. These calculations are controlled by 
input parameters such as launch angle, from which the initial, 
increment and maximum launch angles are used to compute all 
possible ray paths. RAYTRACE "shoots" rays from the 
transmitter and notes where each one "lands" in the vicinity 
of the receiver. The two closest- landing rays are noted and 
the ray which completes the path is calculated by 
interpolation of nearby ray paths. The resulting ray is called 
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a mode. The launch angle is then increased until a second 
mode is obtained. The calculation process terminates when one 
of the following conditions is met: 

* The initial ray does not bounce off the ionosphere but 
rather penetrates it. 

* The launch angle reaches the maximum specified value. 

* The maximum number of modes are reached by the program 
itself . 

The user in ADVANCED PROPHET can define ionospheric 
parameters such as: 

* Critical frequency of the E- region (foE) . 

* foEs if the Sporadic E- layer is to be considered. 

* Critical frequency of the FI region, foFl. 

* Critical frequency of the F2 region, foF2. 

* F2 -layer maximum density altitude, hmF2 . 

* F2- layer semi- thickness , YmF2 . 

When these parameters are specified, the ionosphere will be 
assumed to be uniform throughout the region of propagation. 
If these parameters are not user specified, ADVANCED PROPHET 
uses internal models to calculate the parameters including 
variations along the path due to time of day, solar activity 
and the link geographical region [5] . 

2. ADVANCED PROPHET 4.3 Data Analysis 

After all calculations are completed, two types of 
output displays are generated. The first type of display 
shows the parameters which defined the calculations and the 
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plots of the rays paths, as shown in Figure 18. The second 
display details the results of the calculations in 
alphanumerical fashion, as seen in Figure 19. For the 
comparison process, the rubric, or unique descriptor, 
corresponding to the mode delay (msec) was used to identify 
the different mode path geometries; this delay is the total 
calculated group delay. 

To consider a propagation mode delay prediction valid 
for the comparison, the signal to noise ratio (SNR) should be 
greater than 0 dB. From the alphanumeric output of the 
Raytrace program the SNR can be determined for each of the 
propagation mode delays. The other rubrics of the display 
lists calculated propagation modes, such as number of hops, 
ionospheric bounce region and trajectory losses. For the 
analysis, data obtained from RAYTRACE alphanumeric output were 
utilized. 

C. DATA COMPARISON 

Figures 20 through 25 show experimental and ADVANCED 
PROPHET 4.3 generated mode delay data. Data plotted for a 
given day were assigned "o" and "x" corresponding to PROPHET 
and experimental data respectively. There is a similarity 
between the data, especially in the 2-7 msec time frame. 
PROPHET predicts mode delays that are greater than 15 msec, 
which are not observed in the experiment. In addition, 
PROPHET shows increasing mode delays for the first six hours 
of the day, 0000 through 0600 GMT (1600 through 2200, the 
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Figure 18. 



Ray Paths predicted by ADVANCED PROPHET 4.3 
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*** UNCLASSIFIED *** 










RAYTRACE SYNOPSIS — LAUNCH ANGLES: START- 


.00 END- I 


19.00 INC— 


2.00 


DATE: 10/ 7/92 TIME: 15 


: 00 UT ATMOSPHERIC 


NOISE: NO 1 


BWIDTH: 1 


.240 KHZ 


FREQ: 5.6 


10 cm FLUX: 


137.0 Kp: 2.0 MANMADE NSE: QK 


SNR REQO 


: 12 . OdB 


MONTEREY 


LAT: 36.4 


W LON: 121.6 ANT: 


501 I -OMNI* 


PWR: 


75.00 


SDIEGO 

IONOSPHERE: 


LAT: 32.4 
FOE- 2.5 
HMF2- 344. 


W LON: 117.1 ANT: 500 • .OMNI- 
MHZ FOF1- 3.5 MHZ FOF2- 9.9 
KM YMF2-109 . 8 KM 


RANGE: 


599.3 KM 


NHOP 


1 


1 1 


1 


2 


3 


MODE 


1000000 


1000000 1000000 


3000000 


3300000 


3330000 


ANGLE 


18.65 


23.50 23.90 


45.35 


62.45 


70.80 


DELAY (MSEC) 


2.147 


2.233 2.247 


2.969 


4.509 


6.354 


LOSS ( DB) 


127.29 


147.35 169.46 


120.01 


128.06 


138.42 


GAIN TX/RX 


-4/ -3 


-2/ 0 -1/ 0 


1/ 3 


1/ 4 


1/ 4 


1HZ SHR(DB) 


49.66 


34.04 12.24 


66.52 


60.37 


49.94 


ADJ SNR(DB) 


18.73 


3.11 -18.70 


35.59 


29.43 


19.00 


VIR HTl(KM) 


110. 11 


140.36 143.92 


326.29 


303.14 


302.23 


VIR HT2(KM) 


. 00 


.00 .00 


.00 


302.03 


302.66 


VIR HT3(XM) .00 

RA> 
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. 00 END- 
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:00 UT ATMOSPHERIC 
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2 


MODE 


1000000 
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3300000 


ANGLE 
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DELAY (MSEC) 
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2.364 2.873 
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3.122 


4.438 


LOSS (DB) 


140.94 


135.84 135.38 


176.66. 
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137*59 


GAIN TX/RX 


-4/ *3 


0/2 0/2 


1/ 3 


1/ 3 


1/ 4 


1HZ SNR ( DB) 


35.25 


49.14 50.80 


9.77 


57.66 


50.82 


ADJ SNR(DB) 


4.31 


18.21 19.87 


-21.17 


26.72 
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VIR HTl(KM) 


105.79 
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336.69 


360.24 
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.00 


.00 .00 


.00 


.00 
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.00 
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Figure 19. Alphanumeric output corresponding to 
Raytrace predictions. 
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MODE DELAY vs. REAL TIME (8/28y92) [Prophet (o), Experimental (x)] 
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Figure 20. Comparison between experimental and PROPHET 
Mode Delay predictions for 28 Aug 92. 
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MODE DELAY vs. REAL TIME (10/2^2) [Prophet (o), Experimental (x)] 




Figure 21. Comparison between experimental and PROPHET 
Mode Delay predictions for 2 Oct 92. 

(5.604 - 75 watts) . 



36 



Zulu Hour 




o 

<N 



wn 



O 



3 




N 



Figure 22. Comparison between experimental and PROPHET 
Mode Delay predictions for 3 Oct 92. 
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Figure 23. Comparison between experimental and PROPHET 
Mode Delay predictions for 4 Oct 92. 
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Figure 24. Comparison between experimental and PROPHET 
Mode Delay predictions for 7 Oct 92. 
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MODE DELAY vs. REAL TIME (10/2^2) [Prophet (o), Experimental (x)] 




Figure 25. Comparison between experimental and PROPHET 
Mode Delay predictions for 2 Oct 92. 

(11.004 - 100 watts) . 



40 



Zulu Hour 



local time for San Diego) , which is similar in some aspects to 
experimental results. Single and multihop propagation of up 
to six hops are predicted by PROPHET. The experimental data 
plotted in Figures 12 to 16 shows only three modes. 

In Figure 25 for 2 Oct 92, PROPHET data obtained at 11.004 
MHz with 100 watts is close to the experimental data, where in 
the early hours of the day, propagation is not supported by 
the path. Good correspondence between experimental and 
PROPHET data is noted between 1600 and 2400 GMT. ADVANCED 
PROPHET 4.3 predicts mode delays above 4 msec that are not 
shown in the experimental data. 

For further comparison between experimental and PROPHET 
data, an hour by hour, day by day breakdown is provided in 
Figures 26 through 31. These plots were created by counting 
the times for which the experimental and predicted delays were 
a different by less than 0.5 msec. All hours in which 
propagation for both cases was not available are marked with 
an asterisk. From Figures 26 to 31 it is seen that maximum 
correlation between PROPHET and experimental data occurs in 
the early hours of the day and between 1400 and 1900 GMT. 

The third comparison was calculation of the total 
percentage per hour, over the five day period, in which the 
time delay differences were less than 0.5 msec. In this 
manner a global comparison is performed to more fully 
appreciate how the percentages are distributed throughout the 
day. Table 4 and Figure 32 show this comparison. All of the 
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Figure 26. Mode Delay Approximation Distribution for 
28 Aug 92 
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2. Oct 92 (5.604 MHz - 75 watts) 
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Figure 28. Mode Delay Approximation Distribution for 
3 Oct 92 
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MODE DELAY APPROXIMATION DISTRIBUTION (10/4/92) 




UHUO 



Figure 29. Mode Delay Approximation Distribution for 
4 Oct 92 
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Figure 30. Mode Delay Approximation Distribution for 
7 Oct 92 
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Figure 31. Mode Delay Approximation Distribution for 
2 Oct 92 (11.004 MHz - 100 watts) 
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comparisons and plotted results were made with the aid of 
MATLAB programs . 



TABLE 4. MODE DELAY APPROXIMATION PERCENTAGE DISTRIBUTION 



GMT 


Percentage 


GMT 


Percentage 


0-5:00 


7 . 143 


12 : 00 


1.786 


01:00 


2.679 


19:00 


0.893 


02:00 


4 . 018 


14 : 00 


6.696 


03:00 


5 . 804 


15 : 00 


8.929 


04:00 


3 .571 


16:00 


12.946 


05:00 


1.786 


17:00 


10.268 


06:00 


4 .911 


18:00 


2.678 


07:00 


6.250 


19:00 


0.446 


04 : 00 


3 . 125 


20:00 


0.446 


09 : 00 


4 .911 


21:00 


0.893 


10:00 


2.232 


22 : 00 


2.232 


11 : 00 


2 . 679 


23 : 00 


2.678 
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Figure 32. Mode Delay approximation percentage distribution. 



49 



V. CONCLUSIONS AND RECOMMENDATIONS 



The field measurements indicate that the PROPHET 
propagation predictions are most accurate for the early hours 
of 1400 until 1700 GMT of each day. Data taken on 2 Oct 92 at 
11.004 MHz show that when skywave propagation is supported, 
the experimental and program predictions are almost identical, 
for mode delays of less than 4 msec. As expected, 
experimental and program results are similar also for the part 
of the day when skywave propagation is not supported by the 
ionosphere . 

Field data taken at 5.6 MHz showed good agreement with 
PROPHET propagation predictions. Between 0000 and 0600 GMT, 
PROPHET predicted six mode delay groups (multi -hop 
propagation) with delays that exceed 15 msec. This same trend 
was observed experimentally; however, only three mode delay 
groups were present between 0000 and 0400 GMT. 

It can be concluded that the best correspondence between 
predictions and measurements occurs at 0000 GMT and between 
1400 and 1700 GMT. Reasonable correspondence occurs for the 
hours of sunset and midnight on the Monterey- San Diego path. 

In comparing ADVANCED PROPHET 4.3 predictions with 
experimental field data, it is preferable to input the values 
of the required ionospheric variables corresponding to the 
time the measurements were taken. Since these values were not 
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available for the data used in this research, approximate 
values based on PROPHET'S internally generated parameter set 
were used. 

More data is needed, covering the entire HF band, and will 
be collected as the equipment and software are refined. Once 
the questionable CCIR recommendations for measuring signal 
strength are resolved, SNR can be measured by standard 
procedures and compared with PROPHET predictions. Additional 
ionospheric propagation software such as IONCAP, ICEPAC, ASAPS 
and AMBCOM should be evaluated along with PROPHET. 
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